Nucleon mass: 
from lattice QCD to the chiral limit 
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Previous extrapolations of lattice QCD results for the nucleon mass to the physically relevant 
region of small quark masses, using chiral effective field theory, are extended and expanded in several 
directions. A detailed error analysis is performed. An approach with explicit A(1232) degrees of 
freedom is compared to a calculation with only pion and nucleon degrees of freedom. The role of the 
A(1232) for the low-energy constants of the latter theory is elucidated. The consistency with the 
chiral perturbation theory analysis of pion-nucleon scattering data is examined. It is demonstrated 
that this consistency can indeed be achieved if the A(1232) dominance of the P-wave pion-nucleon 
low-energy constant C3 is accounted for. Introduction of the A (1232) as an explicit propagating 
degree of freedom is not crucial in order to describe the quark mass dependence of the nucleon 
mass, in contrast to the situation with spin observables of the nucleon. The dependence on finite 
lattice volume is shown to yield valuable additional constraints. What emerges is a consistent and 
stable extrapolation scheme for pion masses below 0.6 GeV. 

PACS numbers: 12.38.Gc, 12.39.Fe, 14.20.Dh 



I. INTRODUCTION 

In recent years the nucleon mass has been in the focus 
of steadily improving lattice QCD computations. Techni- 
cal limitations have so far restricted lattice QCD results 
to quark masses larger than five times the physically rel- 
evant masses of the light quarks, m„^d < 10 MeV. Under 
such conditions, the nucleon masses produced on the lat- 
tice are considerably larger [Mm > 1.2 GeV) than the 
physical one. Systematic extrapolations guided by well- 
defined rules of low-energy QCD are necessary in order 
to bridge this gap. Several versions of such extrapola- 
tions have been developed in recent years P, 0, 13 • They 
differ partly in details but agree on the basics, namely on 
the relevance of one-loop chiral pion-nucleon dynamics in 
determining the dependence of Mjv on the pion mass m^. 

The connection between and the u- and d-quark 



masses (for which we take their average. 
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md)/2), is given in leading order by PC AC and the Gell- 
Mann-Oakes-Renner relation, m\f^ — —2mq{qq), with 
the pion decay constant /^r and the chiral condensate (qq) 
both taken at the chiral limit. Accurate lattice QCD re- 
sults for the pion mass as a function of the quark mass 
i^q ^ IH 13 demonstrate that this leading linear relation 
between and iriq remains remarkably stable even at 
large quark masses. Corresponding lattice data are con- 
sistent with one-loop chiral perturbation theory at next- 
to-leading order, up to to^ ~ 0.5 GeV 0,0- Moreover, 
the data continue to display the leading-order (linear) 
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PCAC behavior even beyond this margin, for reasons not 
yet understood in detail. We will therefore accept this as 
a fact when we later compare our theoretical Mjv(TO7r) 
with lattice results. 

The present work updates and extends previous inves- 
tigations 0,01 in several respects: 

• The overall consistency of parameters and low- 
energy constants in the expansion of Mj\j{mT^) with 
those extracted from pion-nucleon scattering, is 
carefully examined. 

• The role of explicit A(1232) degrees of freedom in 
the one-loop nucleon self-energy is studied. 

• A systematic error analysis for the extrapolation 
from the lattice data through the physical point 
down to the chiral limit is performed. 

• A study of finite volume effects in comparison with 
lattice data provides interesting additional con- 
straints. 

Issues of convergence when carrying chiral expansions 
over relatively large ranges of quark masses, will also be 
addressed. 

The following section briefly summarizes the frame- 
work of the present approach, covariant baryon chiral 
perturbation theory (B^PT) using infrared regulariza- 
tion. We omit derivations of basic equations which have 
already been reported elsewhere Q]. Section llTTI describes 
the detailed error analysis of the extrapolation from lat- 
tice data down to the physical region. Section CxI dis- 
cusses connections with the analysis of low-energy pion- 
nucleon scattering. Section introduces the A(1232) as 
an explicit degree of freedom and examines its relevance 
to the chiral extrapolation of the nucleon mass. Section 
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IVII investigates effects induced by varying the finite vol- 
ume of the lattice, and section IVni completes the analysis 
by exploring the effects of avoiding the inclusion of the 
physical nucleon mass as input. Conclusions are drawn 
m section rvml 



II. THE NUCLEON MASS IN CHIRAL 
PERTURBATION THEORY 

In previous work we have investigated the quark 
mass dependence of the nucleon mass using chiral effec- 
tive field theory in the two-flavor, pion-nucleon sector. 
A fit has been performed to lattice data selected accord- 
ing to the largest available lattice volumes and smallest 
accessible lattice spacings, in order to be as close as possi- 
ble to the infinite volume and continuum limits for which 
these xPT calculations are designed. We first present a 
brief summary of this approach. 

In ref. P] a one-loop calculation has been performed, 
taking into account diagrams up to chiral order p"^ in co- 
variant BxPT with infrared regularization j^. The fol- 
lowing result is obtained for MN{m-„) when expanded in 
powers of to^: 
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This formula requires the following input: 



i) the nucleon axial vector coupling constant gA and 
the pion decay constant /tt , both taken in the chiral 
(m,r = 0) limit. In practice we use their physical 
values, gA = 1-267 and = 92.4 MeV, for guid- 
ance and examine variations around these values; 

ii) the nucleon mass in the chiral limit, A/q; 

iii) the low-energy constants ci, C2 and C3. Here ci is 
closely linked to the pion-nucleon sigma term, while 
C2,3 encode information primarily on the A(1232) 
excitation in P-wave pion-nucleon scattering; 

iv) a parameter e^^''(A), combining three different cou- 
plings of the most general 0{p^) BxPT Lagrangian 

(see ref. 0). e.^i\^) represents unresolved short 
distance dynamics. This term compensates the log- 
arithmic dependence on the renormalization scale 
A so as to ensure scale independence of the result. 

The present work builds upon a successful fit of eq. 
to unquenched two-flavor lattice results, referred to as 
"fit 11" in ref. 1]. The input lattice data with improved 



Wilson fermions (points 19, 41, 8 and 23 in ref. and 
table EJ are chosen according to the following criteria: 
small lattice spacing, a < 0.15 fm, and large spatial lat- 
tice size, TTiTri > 5. We have considered only the small- 
est available pion masses, m^r < 0.6 GeV. In order to 
avoid an under-determined fit we fix some parameters. 
In particular, ci is eliminated by substituting the empir- 
ical nucleon mass into eq. JQ) at the physical value of the 
pion mass, C2 is set equal to 3.2 GeV~^ as determined in 
and C3 is fixed at —3.4 GeV~^, according to the NN 
phase shift analysis in ref. ■ Employing the physical 
values /tt = 92.4 MeV, gA = 1.267, only two parameters 
remain free: Mq and e^^''(A) (we chose a renormalization 
scale A = 1 GeV throughout this work) . 

In eq. we have truncated the 0{p^) infrared regu- 
larized expression at in order to avoid further uncon- 
strained counter terms, cf. ref. Q. We note that the trun- 
cation of A/]^-* at 0{m'^) coincides with the 0(j)^) result 
in heavy baryon chiral perturbation theory (HB^PT) , 
see ref. |2| . We stress that in the context of baryon xPT 
with infrared regularization, recoil corrections are part 
of the same order in the corresponding perturbative, di- 
agrammatic expansion. Manifestly covariant and non- 
relativistic methods just lead to a different organization 
of the perturbation theory. We have checked that fits 
based on the HB^PT O(p^) expression give results that 
are compatible with the analysis presented here. Thus 
the term proportional to in eq. i.e. the leading 
recoil correction to the non-relativistic result, plays no 
significant numerical role. 

The application of chiral perturbation theory to inter- 
polations spanning a rather large interval of quark masses 
has always been a point of concern. In the next section 
we examine the degree of convergence as the calculation 
of the nucleon mass evolves order by order in the chiral 
expansion. In preparation of this study we first recall 
further basic results discussed in ref. At chiral order 
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and the one- loop expression at 0{p'^) using infrared reg- 
ularization, expanded in powers of to^, is 
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Here, e^^^(A) absorbs the A-dependence at 0{m^), 
which differs from that of e^'*''(A) in eq. 
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FIG. 1: Nucleon mass as function of m^. Shown is the best 
fit interpolation between lattice results and the physical point 
(star), performed at chiral order using eq. Q (solid curve), 
with input given in column (a) of table |I] The physical point 
is included. The lattice data points in the grey region have 
not been used as input. Also shown are intermediate steps 
at orders and according to eqs. (|5J, J^l- The parameter 
ep'(A) has been fitted to lattice data. 



III. ERROR ANALYSIS 

In this section a systematic error analysis is performed 
in order to test the rehability of the interpolation between 
lattice results and the physical nucleon mass based on 
eq. (HJ. 

We are confronted with uncertainties from two dif- 
ferent sources of errors, namely input uncertainties and 
theoretical uncertainties, which need to be distinguished 
throughout the analysis because they require different 
treatments. Input uncertainties result from statistical 
errors of lattice calculations^ and from limited empirical 
information on low-energy constants. Theoretical uncer- 
tainties come from two sources: First, they arise because 
our fit function derives from a perturbative expansion 
and is therefore subject to corrections of higher order. 
Secondly, we are dealing with an effective field theory, 
which has a limited range of applicability. Our analysis 
is performed in two steps. First, we compare different or- 
ders in the perturbative expansion of Afjv(™ir) focusing 
only on pion masses up to 0.6 GeV. This range is consis- 
tent with the conclusions drawn in ref. In a second 
step we treat input uncertainties. 

Using ci and Mq as given by the M^'' best fit in table 

13 we display, in fig-H m'^^ and examine of eq. 
The pattern of successive steps is certainly promising for 

< 600 MeV. Note that the parameter ef (A) in M^^ 

and e^''''(A) in M^' have different A-dependence. In con- 



trast to ref. m , we do not identify the two parameters at 
a specific scale. Instead e\ (A) is fitted to lattice data, 

yielding ef^(lGeV) ~ 0.48 GeV"^^ a naturally sized 
value for that coupling. A more stringent test of conver- 
gence would have to involve an estimate of corrections 
from Oij)"^). Then, however, the number of poorly con- 
strained parameters would become prohibitively large. In 
the present work a consistency check of the low-energy 
constants required by the fit to the nucleon mass in com- 
parison with the same constants deduced from the analy- 
sis of pion-nucleon scattering data provides a non-trivial 
test that will be performed in section Hvl 

Let us now study the numerical impact of input un- 
certainties for eq. (^). The relevant technical details are 
summarized in appendix ^ The uncertainty in C3 is 
substantial but difficult to quantify. For the moment, we 
ignore uncertainties in C2 and C3 and defer this part of 
the discussion to sections llVI and IVll To start with, 
and are set equal to their values at the physical point. 
Fitting with this setup has led to the results of "fit 11" 
in We now also take into account uncertainties in 
extracting the lattice pion mass, see appendix I A 51 They 
turn out to have little effect. The resulting parameter 
values are listed in column (a) of table |l| The global 
error band for the pion mass dependence of the nucleon 
mass is the "statistical band" displayed in fig.|21 

The values of gA and /^r, to be taken in the chiral 
limit, are expected to differ slightly from the values at 
the physical point. Wc assume to find these values in the 
intervals 

gA = 1.1 ... 1.3, U = 86.2MeV ... 92.4MeV . (4) 

The range for gA is taken from |l3j| . For the lower bound- 
ary of we take the estimate from an analysis of the 
pion mass dependence of in ref. 0| . Varying gA and 

within these intervals extends the error band to what 
is labeled "systematic envelope" in fig. El (see appendix 
The corresponding parameter bounds can be found 
in column (b) of tabled 

In the limit m^r — > 0, the error band remains narrow, 
reflecting our ability to determine Mq quite precisely. Of 
course, this prediction relies heavily on the inclusion of 
the physical point which is located close to the chiral 
limit. For pion masses below 0.6 GeV, i.e. left of the four 
selected lattice points, the band does not bulge much. 
As long as information about the physical point and 
the LECs is included, xPT provides a stable interpolant 
which is well-conditioned to make predictions within its 
range of validity.^ However, the predictive power of our 
analysis in the region to^ > 600 MeV is evidently low. 
The agreement between best-fit curve and lattice data 
up to TO^ « 750 MeV, a scale not small compared to 



At the moment we cannot estimate the intrinsic systematic un- 
certainties in the lattice calculations. 



^ We refer to section 7 for a scenario without inclusion of the phys- 
ical point. 
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TABLE I: Input and output parameters in the fits to large volume lattice data. Lattice pion mass errors have been taken into 
account. The pion-nucleon sigma term is obtained using cr]v = (9Afjv/9m^) with eq. Q. 





(a) statistical error 


(b) systematic envelope 


ei*'(lGeV) 


(GeV-3) 


0.74 ± 0.18 


fitted 


0.27 


... 1.14 


fitted 


Mo 


(GeV) 
(GeV-i) 


0.882 ± 0.003 


fitted 


0.876 


... 0.888 


fitted 


Cl 


-0.93 ± 0.04 


elim. 


-1.04 


... -0.82 


elim. 


9A 




1.267 


fixed 


1.1 


... 1.3 


scanned 


u 


(MeV) 


92.4 


fixed 


86.2 


... 92.4 


scanned 


C2 


(GeV"') 


3.2 


fixed 




3.2 


fixed 


C3 


(GcV-^) 


-3.4 


fixed 




-3.4 


fixed 


xVd.o.f. 


0.13 


0.1268 


... 0.1346 




crjv 


(MeV) 


49 ± 3 


44 


...54 






physical point 
attice data (input) 
. data, m„>0.6 GcV 

fit 

statistical error 
systematic envelope 



0.3 0.4 
[GeV^] 



FIG. 2: Error band at 68% confidence level ("statistical er- 
ror" ) and envelope of bands encoding input parameter uncer- 
tainties ("systematic envelope"). 



the chiral symmetry breaking scale A;^ w ^T^fir, has low 
statistical significance.'^ 

Note that while the individual fit parameters in tabled 
exhibit appreciable sensitivity to the input choice for 
and /tt, the overall shape of the fit curve shows hardly 
any dependence on gA, /tt in the range of interest for 711-^, 
such that the "systematic envelope" in fig.[21for the global 
error band is barely distinguishable from the band found 
for fixed gA and 

In ref. [Tsl l an attempt was made to determine the im- 
pact of higher-order effects by plotting a band for 
selecting the chiral limit values of gA and /^r from error 
intervals while all other parameters were kept fixed. How- 
ever, when studying higher-order effects one must sub- 
stitute the pion mass dependent functions gAi^n-n) and 
/^(m.n-), constrained in such a way that they are consis- 
tent with lattice QCD results over the whole range of m^. 
This has not been done and therefore the significance of 



^ Complementary results for large quark masses are obtained using 
diflerent methods, such as the Adelaide approach 3j and the 
chiral quark soliton model 1341 . 



the analysis [13 remains doubtful. 

Our conclusion is instead that the interpolation based 
on 0{p'^) chiral perturbation theory is sufficiently stable 
for < 0.6 GeV. 



IV. CONSISTENCY WITH PION-NUCLEON 
SCATTERING 

In order to test whether our extrapolations from lat- 
tice results are physically meaningful, we compare them 
with pion-nucleon observables for which chiral expan- 
sions involve the same parameters but in different com- 
binations. One such observable is the isospin-even t:N 
S-wave scattering amplitude at threshold, T+(m7r) = 
47r (1 -t- m-^/MN) a"*" where a"*" is the corresponding scat- 
tering length. Empirically one finds the anomalously 
small value T+ = (-0.12 ± 0.11) GeV"^ Chiral symme- 
try implies that the leading term of at chiral order p 
vanishes. In HBxPT the contributions at order p^ and 
p^ read 



2mi 
P. 
0{p') 



C2 + C3 - 2ci 



A. 

8M0 



6471/4 



(5) 



The term of order p'^, supposedly small, involves a series 
of additional low-energy constants which arc not well de- 
termined and therefore presently of no practical use for 
a detailed estimate. 

Another suitable observable is the isospin-even, spin- 
averaged P-wave ttN threshold amplitude, P^{mTr) = 
47r(l -I- m^/Mjv)(4a33 + 2a3i -t- 2ai3 + an), given in terms 
of the P-wave scattering volumes 02/, 2J in channels with 
spin/isospin J, I = 1/2 or 3/2. Empirically, one finds 
Pi+ = (1044 ± 38) GeV-^Jn (l^. In HBxPT, the terms 
up to chiral order p'^ are 




2 

II 
0{p') 



C2 



C3 



g^m,^ 



g\m^ 



Ap^M^ 12^/4 



77 

32 



(6) 



where one should note that Pi enters in the scattering 
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amplitude with an extra factor involving pion momenta, 
q-q', of order in the chiral counting. 

Further independent information can be drawn from 
the pion-nuclcon sigma term, 



(7 A/- — (A^|m„ uu + rud dd\N) — m, 



DM 



N 



dr 



(7) 



Making use of the Gell-Mann-Oakes-Renner relation, 
one obtains: 



2 dMN . 2 

omt 



(8) 



Ref. [13 provides us with the empirical value aN = (45 ± 
8) MeV. 

The requirement that the empirical T+, Pj^and ctat are 
simultaneously reproduced imphes C3 ~ —(5 ± 1) GcV^^ 
if one uses the 0{p^) HBxPT formulae eqs. (O, (jHl at 
the physical pion mass. It might then appear that C3 
determined from ttN scattering is incompatible with the 
value C3 = — 3.4GeV~^ consistent with NN scattering 
and used as input previously. At tree level and in the 
non-relativistic limit, the spin-isospin averaged P-wave 
scattering volume is well known [23, |^ to be propor- 
tional to (A^ — ti>^)^^, where the pion energy uj equals 
TOtt at threshold. The fact that the delta-nucleon mass 
difference, A = Ma — M^, is barely twice the physical 
pion mass, m^^^^, is at the origin of the strong energy de- 
pendence of the P-wave TriV amplitude. A determination 
of C3, first by comparison with ttN threshold data and 
secondly by examining its role in peripheral NN phase 
shifts, will therefore lead to different conclusions concern- 
ing C3, the NN situation being effectively closer to the 
limit of a static pion field with a; = |22j | . This apparent 
discrepancy just reflects the in-effectiveness of HB^PT 
to deal with the well-known strong energy dependence 
of the ttN amplitude. The important effects of this en- 
ergy dependence are then "hidden" by absorbing a large 
correction of order m^/A^ into C3 when using the 0{p^) 
result of eqs. jSJ, © for its determination. Including this 
leading correction the estimate of C3 is reduced by a fac- 
tor of 1 — (m7r/A)^ « 3/4, which accounts for much of 
the difference between the two cases. The large value of 
C3 « —5 GeV~^ is therefore an artifact of the HB^PT 
expansion truncated at 0{p^). 

Matching the tree-level ttN amplitude with explicit 
A(1232) and the ttN ttN graph calculated from 
the second-order pion-nucleon effective Lagrangian, the 
A(1232) contribution to C3 reads 



44 
9A 



(9) 



where ca = SttJVA /7r/(2Mo) in terms of the ttA^A cou- 
pling constant g-r^NA, and all constants (/7r,A,Mo) are 
understood to be taken in the chiral limit. With the 
frequently used empirical coupling ca — 1-5, see section 
IV Bl we have ~ —3.4 GeV~^, which agrees with C3 ex- 
tracted from the NN scattering analysis. We interpret 



this value as being representative for the C3 to be used in 
the extrapolation of lattice results for Afjv- 

We point out here that the set of values Mq — 



0.883 GeV, ci 



-0.89GeV-\ C2 = 2.98GeV~\ C3 



-3.55GeV-i, e^'*^(lGeV) = 0.46 GeV^, gA = 1.1 and 
/tt = 92.4 MeV represents an optimal fit to lattice data, 
gives a curve right in the center of our statistical error 
band fig.El and satisfies the empirical constraints for cr^Vj 
and provided that C3 in eqs. |3|), is replaced 



by (4/3)c3 



We note that the factor 4/3 is to be 



understood as representing the substantial contribution 
at C'(p^) in eqs. lO, ® arising from the A(1232) propa- 
gator structure. Our intermediate conclusion is that the 
parameters required by the "best fit" to lattice results 
for Mat are consistent with those from the analysis of ttN 
and NN low-energy data once C3 is corrected for artifacts 
of the Heavy Baryon truncation. These considerations 
will be strengthened as we now further examine the role 
of the A(1232) and investigate the systematics of finite 
volume effects. 



V. EXPLICIT A (1232) DEGREES OF FREEDOM 

An effective field theory which includes only pion and 
nucleon as explicit degrees of freedom encodes implicitly 
contributions from the A(1232) resonance through low- 
energy constants, e.g. C3, but its interpretation requires 
caution, as elaborated in the previous section. Work- 
ing at limited perturbative order, "freezing" the A(1232) 
and relegating its effects to higher order terms can lead 
to a rather in-effective chiral expansion, given that the 
delta-nucleon mass difference is only about a quarter of 
the chiral symmetry breaking scale, A;^ w '^T^fir- It has 
been shown that including the A(1232) as an explicit 
degree of freedom in spin-dependent quantities like the 
magnetic moments or the axial coupling of the nucleon 
promotes important quark-mass dependent contributions 
to low orders in the perturbative calculation, leading 
to well-behaved chiral extrapolation functions, see e.g. 
refs. O nil' For the case of the nucleon mass explicit 
A(1232) treatment turns out to be less crucial. However, 
it helps us in clarifying the role of the p^-effects in eq. , 
which are dominated by the couplings c^. 



A. Formalism 

We work in the Lorentz covariant formulation of the so- 
called small scale expansion (SSE) introduced in ref. [23 |. 
In this scheme the delta-nucleon mass difference in the 
SU{2) chiral hmit, A = Ma - Mq, is included in the 
power-counting, together with the pion mass and soft 
external momenta, as a small scale generically labeled e 

m 

The main topic of our discussion with explicit delta 
degrees of freedom is to clarify the success of the C'(p^) 
calculation without A(1232). The leading effect due to 



6 




FIG. 3: Leading-one- loop diagram for the nucleon self-energy 
with an intermediate A (1232). 



explicit delta appears at order e"^ in SSE and starts to 
contribute at to^, like the 0{jp^) graphs. Therefore we 
have concentrated on the 0{e^) leading one-loop contri- 
bution to the nucleon self-energy involving the propaga- 
tion of the A(1232), fig.El'* In Section rvO we will com- 
pare our results with the expression in ref. 26] in 
order to estimate the importance of higher-order effects. 

The €^ calculation requires the leading chiral pion- 
nucleon-delta effective Lagrangian 



CttNA = 



/it 



h.C. 



(10) 



where 4'^ is the Rarita-Schwinger field representing the 
A(1232) and is the nucleon field. We employ the fol- 
lowing propagator for the free spin-3/2 isospin-3/2 field 



i + Af A 



p3/2 



(11) 



mass of the nucleon, cf. l2£ 



5M^ = 



id -2) 



4f^'{d-l)MoMl 



- Ml - m2)(M2 + 2AfAAfo + Ml - m2)2 



A/ 2 M 2 



2A/AAfo + 4m^ 



d-l 



{Ml-ml){2MAMo 



Ml 



(12) 



where d = 4 is set after removal of the singularities from 
I A and A^r- The resulting expression for the nucleon 
mass at order is then of the form 



Mn = Afo - 4 cimi + SM^"" 
-\- counter terms , 



SM 



N 



(13) 



where 6MJ^^ is the nucleon mass-shift at order in the 
scheme with ttN degrees of freedom only, cf. eq. We 
note that we recover the non-relativistic 0{e^) SSE result 
for the nucleon mass^ by just keeping the leading term 
in the I/Mq expansion in eq. (|13|l with eq. (|12|l inserted. 

In a first approach we truncate at order 1 /Mq the result 
of eq. (221. In addition to the term — 4 ca ^at^'at, 
the SSE-scheme provides two more counter terms at this 
order'' 



where P^l^ and ^g'^j are the spin- and isospin-3/2 projec- 
tion operators, respectively. Using the propagator pi|l . 
the computation of the one-loop A(1232) contribution to 
the nucleon self-energy in infrared regularization requires 
the standard loop integrals ^ 



,(3) 
-ttN 



[Bi A^ + 4B2ml A] * 



(14) 



which are sufficient for renormalization. The finite parts 
of the renormalized couplings Bi 2 are then set such that 
the chiral expansion of MN{mTr) at order e'^ begins with 



lA^ = 



(271^ {ml -k^- ie)[Ml - (p - fc)2 - ie] ' 
d'^k 1 



(27r)'^ {ml -k^ - ie) 



M. 



N 



Afn 



Aciml 



0{ml) . 



Imposing this condition, one obtains 



(15) 



Here d denotes the space-time dimension. Both previous 
integrals diverge as d — > 4. 

The 0{£^) graph in fig.O leads to a correction to the 



Bl{X) = 

= 



9(47r^)2 
(47rA)2 





A 

" Aio 


)- 


A 




A 


Mo ^ 


(- 


Aio 



Mo J A 



2A' 



(16) 



* For a leading one-loop calculation in a different framework see 
ref. [35 1 . 

^ See ref. 9] for details. In distinction from this reference, we do 
not identify the regularization scale with Mq, which is a param- 
eter in our approach. 



® See, e.g., the calculation in ref. l36l . which in addition contains 

an "effective" coupling, like e'^^ in section ITTl 
^ The constants Bi and B2 used here are denoted B:j2 and B23, 

respectively, in ref. l26l . 
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It follows that 

Mn ^ Mo-^ciml 



3mz 



327r/2 64 7r2/2Mo 

2 I in„2 



i9A + 3ci) 



(3gi + 104) In ^ + 4 e^{X)mt 
327r^/^ Mo A 



1 



(2A3 - 3m2 A) ( 1 
A 



r2 



67r2/2 



A 



r2 



37r2/2 



1 



2Afn 



2Mo 

(A^ - )3/2 



In — 
2A 



X In ( A + /A. _ 1 
\ TO^ V ml 



(17) 



The last formula is valid for < A. Since in our numer- 
ical analysis we take as input lattice data at pion masses 
larger than the physical one, we need also the analytic 
continuation of the expressions above to TOtt > A. This 
is done with the replacement 




In the next section we will also analyse the non- 
truncated O(e^) result of eq. (|13|l with eq. (|12|l . keeping 
the full tower of 1/Mq recoil corrections. We introduce 
the same counter term structures as for the truncated ex- 
pression and accept uncompensated regularization scale 



dependence < 



and ^ in this case. The finite parts 



of the renormalized Bi are again determined by imposing 
eq. (fTB|l . 

B. Numerical results 

Using as input our selected set of lattice data dis- 
cussed in section Ull we now analyze the 0{e^) expres- 
sions for MNim.^) introduced in the previous section: 
the (!?(l/Afo)-truncated eq. H17II and the "full" expression 
with its residual scale dependence. We fix the renormal- 
ization scale again at A = 1 GcV. The mass splitting 
A is identified with its value at the physical pion mass. 
The axial vector nucleon-delta coupling ca is treated as 
input quantity, determined from the decay width of the 
A(1232) as follows. In the rest frame of the decaying 
delta we have the relativistic expression: 



AMI 



(18) 



where q is the pion or nuclcon momentum in that frame. 
At the delta mass and width corresponding to its pole 
position W\, Ma - iTA/2 = (1210 - i 50) MeV, one 



1.8 

1.6 

> 1.4 
u 
O 

^ 1.2 
1 

0.8 



physical point 
lattice data (input) 
1. data,m;r>0.6GeV 
Stat, err., w/o expl. A 

- - fit delta I 

— fit delta II 



0.1 



0.2 



0.3 0.4 
m/ [GeV^] 



0.5 



0.6 



FIG. 4: Best fit curves based on the formula at order e in 
SSE. The short-dashed curve refers to "fit delta I" , while the 
long-dashed curve corresponds to "fit delta 11" in table ^ 
For comparison, we plot the 68% statistical error band of the 
0{p'^) BxPT result shown in fig.H 



finds Ca — 1.5. We fit again the three remaining unknown 
parameters: Mq, ci and eA(l GeV). 

In table ini "fit deha I" refers to eq. ijlTl, while "fit 
delta 11" refers to the relativistic untruncated expression. 
We fit to the four data with m^r < 600 MeV including 
the physical point, see fig.0] In "fit delta 11" the output 
parameters Mq and ci remain stable (within error bars) if 
we vary the input scale A between 0.5 and 1.5 GeV. Thus 
the residual scale dependence is indeed under control. 
In light of the discrepancies between the two fits, the 
truncation of eq. (|17|l at order 1 /Mq neglects terms which 
seem to play a significant role. Furthermore, we have 
checked that the 68% error bands associated with "fit 
delta I" and "fit delta 11" do not overlap for < 300 
MeV, except at the physical point used as a constraint. 
Therefore we can conclude that the discrepancy between 
the truncated and untruncated expressions does not have 
a statistical origin. This situation is different in the 0{p'^) 
BxPT calculation discussed in ref. 0, where truncating 
at 1 /Mo represents a very good approximation to the full 
result. 

In fig.Qlwe plot the curves corresponding to "fit delta 
I" and "fit delta 11", together with the 68% error band 
of the 0{p*) BxPT result in section lTTll previouslv drawn 
in fig.m At that confidence level, BxPT at Oip^) and 
0{e^) covariant SSE are compatible for the whole range 
of pion masses under study. At the present level of ac- 
curacy, treating the A (1232) as an explicit, propagat- 
ing field is therefore not essential for a satisfactory de- 
scription of the quark mass dependence of the nuclcon 
mass. The virtual A(1232), being far off-shell, has a short 
propagation length. Therefore an equally successful in- 
terpolating function can be obtained by "freezing" the 
delta effects into low-energy constants, working at fourth 
chiral order, with a value for the coupling C3 governed 
by delta-dominance, as discussed in section IIVI Note 
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TABLE II: Fit results for M]v(m.^) in covariant SSE, at leading-one- loop order. "Fit delta I" refers to the formula trun- 
cated at order 1/Mo, "fit delta 11" to the "full" expression with residual scale dependence. The input parameters are 
ca = 1.5, Pa = 1.267, /,r = 92.4 MeV and A = 271.1 MeV. The renormalized effective coupling ba entering at 0{m^) 
has a different regularization scale dependence in the two cases. 





Mo (GeV) 


ci (GeV') 


eA(lGeV) (GeV~^) 


xVd.o.f. 


fit delta I 


0.894 ± 0.004 


-0.76 ±0.05 


4.5 ±0.1 


0.19 


fit delta II 


0.873 ± 0.004 


-1.08 ±0.05 


2.8 ±0.2 


0.43 



that the ttA loop integral in fig.Ol involves a denomina- 
tor {y/q"^ + ml + A)~-^, characteristic of the crossed (u- 
channel) delta pole in the pion-nucleon scattering ampli- 
tude, for which the sensitivity to A turns out to be less 
significant. 

These conclusions hold for any choice of the input pa- 
rameters gA, /it, Ca and A within phenomenologically 
meaningful limits. In order to check the sensitivity with 
respect to variations of the delta-nucleon mass difference 
we have performed fits for A = 293 MeV and 330 MeV. 
The former corresponds to the 90" ttN phase-shift in the 
spin- 3/ 2 isospin-3/2 channel, while the latter is moti- 
vated as a chiral limit value in ref. l2a. 



be we have made a more conservative estimate, varying 
them between — 3...3GeV^^. Scanning those ranges, 
we have performed fits with three free parameters Mq, 
ci and ei(A = 1 GeV). We have fixed gA = 1.267, /^r = 
92.4MeV, A = 271.1 MeV and ca = 1.5. The resulting 
band of best-fit curves is shown in fig.El In output we get: 
Mo = 0.85.. .0.93 GeV, ci = -1.32 ... -0.25 GeV^ and 
ei(lGeV) = -7.3...1.9GeV~^ In view of these results, 
our e'^ results are potentially affected by important sys- 
tematic effects from higher orders. At present, however, 
lack of information on SSE couplings does not permit us 
to perform the analysis discussed in previous chapters at 



C. ©(e^) vs. Oie" 



VI. VOLUME DEPENDENCE 




0.2 0.3 
[GeV^] 



0.4 0.5 



A. Implementing finite lattice sizes 

So far we had to restrict ourselves to the largest avail- 
able lattice sizes. However, xPT is also able to describe 
deviations from the infinite- volume limit due to the finite 
spatial extent of the simulation volume Com- 
bining the L- and m^-dependcncc in our analysis of the 
nucleon mass enables us to enlarge the input data base 
for our analysis. 



FIG. 5: Comparison between our O(e^) "fit delta 11" (dashed 
curve) and a band of 0{e^) fits using natural size assumptions 
for higher-order couplings. 



Eq. (20) of ref. ,2^] provides the contribution to 
Mjv(?Ti,r)-^ This expression involves the poorly known 
SSE couplings C2, C3, 63, bg and ei. Note that C2, 
C3 and ei have different numerical values in SSE and 
BxPT. Motivated by delta-dominance, we assume the 
range —1 ... 1 GeV~^ for C2 and C3 within SSE. For 63 and 



The reader should be aware of misprints in eqs. (19), (25) and 
(26) of that paper. 



(a) 



FIG. 6: One-loop graphs of O(p^) (a) and 0{p'^) (b, c) con- 
tributing to the nucleon self-energy in BxPT (without ex- 
plicit A (1232) degrees of freedom). The solid dot (diamond) 
denotes a vertex from first (second) order nN Lagrangian. 



With pion loop momenta restricted to the discrete 
values permitted by periodic boundary conditions, the 
0(p4) BxPT result for the mass difference SM^ between 
finite and infinite volume has been published in ref. @] . 
At order one obtains from diagram (a) of fig.El using 
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infrared regularization 

(5^3) Ma, 



167r2/2 





/"OO 


f dx + 




Jo 





E 

neZ3\{0} 



i^o (^i|n|^M2x2 + m2(l-a;)^ , (19) 

where the integral from 1 to infinity is taken in accor- 
dance to the prescription in rcf. 9J. 

The 0{p'^) contribution stems solely from the tadpole 
graph diagram (b) of fig.El and reads 



4^ 



C2- 



nGZ3\{0} 

K2{L\n\mT, 



C3j 



L\n\mTr 



(20) 



Here Kq, Ki and are modified Bessel functions. 
Eqs. (fli^ . (pn|l have been worked out in the so called "p- 
regime", namely for rriTrL 3> 1. In ref. Q the parameters 
of eq. iQJ have been fitted for large volumes, L ~ 2 fm. 
The resulting predictions for the finite volume effects at 
0{p'^) showed a surprisingly good agreement with dy- 
namical improved Wilson data of the QCDSF, UKQCD 
and JLQCD collaborations, even down to L « 1 fm. Here 
we use eqs. (|19|l . (I2U|I for a global fit of lattice data in dif- 
ferent simulation volumes. 

We note that no new parameters enter in the finite- 
volume corrections SM^, and that ci, C2 and C3 appear in 
combinations different from those in the infinite volume 
formula Thus the m,r- and L-dependence comple- 
ment each other in constraining these parameters. The 
nucleon mass in the finite volume reads 

MM{m^,L) = MN{m^) 

+ <5(3)Mjv(m^, L) + S^'^^MN{m^,L) , (21) 

where Mjs/^mT^) is the order result of eq. corre- 
sponding to the limit L — + 00. For the numerical evalu- 
ation of 5Mjq{mT^,L), it is precise enough to terminate 
the infinite sums at \n\ = 8 and |n| = 6 in eqs. 119|) and 
(|20(l , respectively. The integral in eq. (|19() is computed 
numerically. 

Appendix ^ lists the lattice data used for this pur- 
pose and explains our selection criteria. The mass 
m eq. (Ell) is the pion mass in the infinite volume. For 
each set of lattice simulations, we have identified the in- 
finite volume pion masses with to^ at the largest volume 
(L > 2 fm) . Eq. H21|) is used to determine low-energy pa- 
rameters by fitting to finite volume lattice data following 
the same statistical strategy as described in section iHll 
We can now afford to release C3 and determine it from the 
fit. Also, we accommodate an uncertainty about C2 in a 
range from 3.1 GeV~^ to 3.3GeV~^, which encompasses 
results from several HB^PT fits at 0{p^^ to experimen- 
tal ttTV scattering data, see table 4 in (29|. The output 
parameters listed in table IIIll are then used to draw 
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physical point 




lattice data, large L 




lattice data, small L 




M,v(m„), L=co 




statistical error 


■ 


systematic envelope 



L = 1 .02 ftn 
1.05 ^LQ2 

LIO 

1.20 ^ T 

1.30 

Tj.58 




0.2 0.3 
m„'^ [GeV^] 

FIG. 7: Finite volume lattice data and inferred nucleon mass 
MNirriTy, L). All lattice points shown here are used as input 
for the fit, cf. table in appendix|n| The solid curve, the 68% 
statistical error band and the systematic envelope project the 
full information, including the finite volume data, onto the 
infinite volume limit, eq. Q. The calculated finite volume 
dependence is shown by the dashed curves. 




1.2 1.4 1.6 1.8 
L[fm] 

FIG. 8: L dependence of the fit function and error bands for 
fixed pion masses. The horizontal dotted lines are the infinite 
volume limit of the fit function. The dashed lines represent 
the ©(p^) contribution to the finite volume correction. 



• best-fit curve and 68% confidence level error bands 
based on the infinite- volume expression eq. JQ), see 

fig-IZI 

• i-dependence at fixed based on eq. (|21|l . see 

fig-El 

Remarkably, C3 comes out low in magnitude, compatible 
with NN scattering results [13 , the outcome of ref. jj] , 
the finite volume study of ref. fg and delta-dominance. In 
the previous infinite- volume analysis, it could be argued 
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1.4 1.6 
L [fm] 

FIG. 9: The box length L dependence of the nucleon mass 
in relativistic SSE 0{e^) compared to the BxPT C'(p*') result 
and lattice data at fixed pion mass = 545 MeV. 
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FIG. 10: Comparison of 4sE^Jv('^'r, in eq. 12211 and 
in eq. I|2()|l . 



that C3 and the counter term parameter e^''^(A) can be 

played against each other, given that e^^^ remains essen- 
tially unconstrained. The L-dependence now provides a 
further condition which definitively points to the "small" 
C3 discussed previously in section HVI 



B. Finite size effects with explicit A(1232) 

In order to check the influence of the A(1232) reso- 
nance also in the finite volume case, we present a calcu- 
lation of finite size effects in SSE at O(e^) with infrared 
regularization. 

The basic ingredients have already been given in sec- 
tion The contribution from the diagram, fig.|3l to 
the difference SAIn between finite and infinite volume 
becomes ISOl: 



c(3) 



1-x 





/•oo 


[ dx-i 




Jo 





E 

nGZ3\{0} 



Ma 
Mo 



y{x)Ko [L\n\^/y{x] 



Llrtl 



K, {L\n\^/y{x) 



(22) 



where y{x) 
we refer to 



l)x- 



-Mlx' 



- (1 — x) , for details 



The finite volume mass-shift at order e'^ 



is then given as: 



(23) 



It has to be added to the infinite volume result of eq. (|13|l . 
Here we do not perform any truncation in the expansion 
in powers of I/Mq. Eq. (|23|l is used with the parameters 
from "fit delta 11" to predict the finite volume effects of 
the nucleon mass with explicit A(1232). 

In fig. Owe show a comparison between SSE at 0{e^) 
and BxPT at 0{p'^) with the parameters listed in tabled 
Both curves are compatible with lattice data. As in the 
infinite volume case, at the present level of accuracy, the 



inclusion of explicit A degrees of freedom is not essen- 
tial to reproduce the trend shown by the lattice data, in 
the sense that effects of the A can be absorbed in low- 
energy constants at order p'*. In order to quantify the 
contribution of the A(1232) to the tadpole diagram in 

BxPT, we compare 4se^^a^ eq. (EH and S^'^^'Mn in 
eq. (1201) for m^r = 545 MeV, see fig-El The overall re- 
sult gives additional support to our arguments in section 
IIVI emphasizing the delta-dominance in the low-energy 
constants entering the tadpole at 0{p^). 

In ref. 'si'l formulae for the finite size effects of the 
nucleon mass have been presented in the framework of 
non-relativistic chiral effective field theory with explicit 
A(1232) degrees of freedom. We have checked that for 
pion masses larger than 300 MeV and L < 1.5 fm the 
corresponding mass difference between finite and infinite 
volume lies systematically below the covariant SSE result 
of eq. 1221), suggesting that recoil corrections could play 
an important role in this regime |30| . 



VII. SYNTHESIS 

With two independent variables at hand, the pion mass 
TO^ and the lattice size L, one can make use of an en- 
larged data base to extrapolate down to small quark 
masses. Figure ITTI shows the results of such a study in 
the framework without explicit A(1232). Now the phys- 
ical point is not included as a constraint. Instead, the 
low-energy constants Ci have been fixed as input. Two 
scenarios have been chosen for comparison, with two dif- 
ferent sets of Ci, as compiled in table Hvl In both cases 
the finite- volume data with L < 2 fm have been included 
in the fit procedure, in addition to the large-L data shown 
in fig. llll The spin-isospin averaged S- and P-wave ttN 
threshold amplitudes are evaluated in parallel, using the 
rule that C3 is saturated by A(1232) contributions, as 
elaborated in section 4. The set denoted "scenario 2" 
evidently meets the constraints imposed by pion-nucleon 
scattering data quite well. The fits show that Mq and e^^^'' 
are well-determined, and the narrow uncertainty band 
permits making precise statements at the physical pion 
mass. Note that an exact knowledge of the Ci is cru- 
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TABLE III: Input and output parameters for the finite volume fit. 



ei'*^(lGeV) 

C3 

Mo 

Cl 

U 

C2 

xVd.o.f. 
aN 



(a) statistical error 



(b) systematic envelope 



(GeV-3) 
(GeV-i) 
(GeV) 
(GeV-i) 

(MeV) 
(GeV-^) 

(MeV) 



0.60 ± 0.09 
-2.9 ± 0.6 
0.884 ± 0.006 
-0.88 ± 0.09 
1.267 
92.4 
3.2 
0.75 
47.6 ± 4.7 



fitted 
fitted 
fitted 
elim. 
fixed 
fixed 
fixed 



n 99 




lib Leu 


-4.3 . 


. -1.4 


fitted 


0.873 . 


. 0.898 


fitted 


-1.03 . 


. -0.69 


elim. 


1.10 . 


. 1.30 


scanned 


86.2 . 


. 92.4 


scanned 


3.1 . 


. 3.3 


scanned 


0.69 . 


. 0.82 




36 . 


. 57 





1.6 



> 

a 1.2 



•ir physical point 
— ^ lattice data, large L 
BBI fit & band, scenario 1 
fit & band, scenario 2 




0.2 0.3 
m^^ [GeV^] 

FIG. 11: Chiral extrapolation of finite volume lattice data 
down to the physical pion mass, based on the 0{p'^) expres- 
sion in eq. (Q. The relevant parameters have been fixed as in 
table HVI Shown are the resulting 68% local statistical error 
bands for the nucleon mass extrapolated to infinite volume. 
The band width can be directly interpreted as the error of a 
prediction at each given pion mass, see appendix lA 41 Only 
those lattice data points corresponding to large volumes with 
L > 2 fm are displayed. 



cial to achieve this resuh. Therefore, in order to release 
those parameters and keep precision, several observables 
involving those low-energy constants should be simulta- 
neously analyzed - including finite size effects. Such a 
combined approach promises to be the winning strategy 
in performing reliable chiral extrapolations. 



VIII. CONCLUSIONS 

In this work chiral extrapolations of lattice QCD re- 
sults for the nucleon mass have been extended in several 
directions, with the following conclusions to be drawn: 

i) A detailed error analysis demonstrates the statistical 
significance and stability of these extrapolations for pion 
masses below 0.6 GeV. 

ii) For the results to be physically meaningful, it is im- 
portant to verify that the low-energy constants of the chi- 
ral pion-nucleon effective Lagrangian which control this 
extrapolation are consistent with those extracted from 



TriV and NN observables. This is indeed demonstrated 
to be the case, provided that the A(1232) dominance of 
the P-wave TriV parameter C3 is properly accounted for. 

iii) Unlike the situation with spin observables such as 
the nucleon axial coupling constant qa, at the present 
level of accuracy, it is not crucial to introduce the 
A(1232) as an explicit, propagating degree of freedom 
when dealing with the mass of the nucleon. The off-shell 
propagation length of the A in the ttA loop correction to 
the nucleon mass is sufficiently short that its effects can 
be absorbed in low-energy constants, with the caution 
exercised as stated in ii). 

iv) A highly useful additional source of information 
from lattice QCD results, apart from their quark mass (or 
m^n-) dependence, is their variation with the finite lattice 
size L. The observed systematics of this variation sup- 
ports the underlying chiral dynamics framework, in terms 
of pion-nucleon degrees of freedom, for the extrapolation 
to small quark masses. This systematics also reduces 
ambiguities in constraining input low-energy constants. 
In particular, the P-wave ttN parameter C3 preferred by 
the chiral extrapolation of the nucleon mass turns out to 
be compatible with the one deduced from NN scattering 
phases shifts. The larger value found from fits based on 
0{p'^) Heavy Baryon xPT to threshold pion-nucleon scat- 
tering data can be explained in terms of pion-nucleon- 
delta dynamics in the P-wave channel. 

In summary, chiral effective field theory extrapolations 
of the nucleoli mass from lattice QCD to the physically 
relevant region of small quark masses and further on to 
the chiral limit are beginning to reach a high degree of 
consistency with other independent low-energy observ- 
ables. In order to further improve the accuracy of such in- 
vestigations, it is an important task for the future (apart 
from expanding the lattice QCD data base to smaller 
pion masses) to perform simultaneous systematic extrap- 
olations of several observables using a single consistent 
set of low-energy constants in the underlying effective 
Lagrangian. 
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TABLE IV: Input parameters and resulting output for the nucleon mass extrapolation based on finite volume lattice data. 
Also listed are the thereshold pion- nucleon amplitudes and . Consistency with empirical values for these observables is 
achieved by substituting cs (4/3)c3 in their O(p^) HBxPT expressions, as discussed in section HVl 





scenario 1 


scenario 2 






(GeV~^) 


0.38 ± 0.04 


0.62 ± 0.04 


fitted 


Mo 


(GeV) 


0.837 ± 0.017 


0.924 ± 0.017 


fitted 


Cl 


(GeV-^) 


-1.0 


-0.95 


fixed 


C2 


(GeV-i) 


3.2 


3.0 


fixed 


C3 


(GeV-i) 


-3.4 


-3.65 


fixed 






1.2 


1.2 


fixed 




(MeV) 


90 


90 


fixed 


UN 


(MeV) 


55.42 ± 0.13 


49.73 ± 0.14 






(GeV) 


0.899 ± 0.017 


0.981 ± 0.017 




T+ (C3 ^ (4/3)c3) 


(GeV"-') 


2.984 ± 0.021 


0.101 ± 0.017 




P+ (C3 (4/3)c3) 


(GeV-3) 


949.4 ± 3.0 


1010.5 ± 2.3 





APPENDIX A: DETAILS OF ERROR ANALYSIS where w'^^ {t) is the cumulative distribution for n 

degrees of freedom. 

We have checked that, for our purposes, is ap- 
proximately quadratic in p around the minimum, i.e. 

^2j-pcst ^p-j ^ ^2j-pcst-j _|_ i^p^ ^^^^ each pion 

mass X, we determine the minimal and maximal val- 
ues of the fit function, /min(2;) and fmaxix), which can 
be found in the confidence region i?cL = + Sp : 

^SpkAk^iSpi < XcL.nl • Shading the region between 
/min(a;) and /max(a^) results in a global "statistical error 
band" such as shown in fig.|21 It sweeps over all inter- 
polating functions allowed by statistics within the confi- 
dence level. 

The error matrix E'^^'^ and the gradient v{x) of the 
interpolating function f{x) are defined as 



1. Confidence regions and global error band 

The information provided by the N lattice data points 
is of the general form 

"Mat at m-^ = Xj ± Axj has values yj ± Aj/j" 
for j = l,...,iV. 

For simplicity, let us first ignore the errors Axj. We 
assume the errors in the Uj to be uncorrelated and dis- 
tributed normally, one standard deviation being Ay^ . We 
denote our fit function f{x;p,u), with parameters col- 
lected in p — (pi,...,p„) and u = {ui, ...,Um)- The pa- 
rameters p are those which we will estimate from lattice 
data. For the other parameters u, most prominently gA 
and /^, we take estimates from literature. During the 
statistical analysis of uncertainties in p, we keep u fixed. 
Therefore, we will omit u in our notation. 

Having agreed on a confidence level CL G (0,1), 
we are looking for a confidence region Rcl fulfilling 
P e Rcl] = CL , i.e. the probabihty to find the 
true parameters p^^'^^ in the confidence region is CL. 
Throughout this work we choose CL ~ 68% which cor- 
responds to one standard deviation for Gaussian dis- 
tributed errors. 

As described in 0, H^, , the x^ method offers a way 
to construct such regions in good approximation. Using 

X^iy\p) = iVj - fixj;p)f /Ay/ we choose 

Rcl ■■={?■■ X^iy\p) < Xopt(2/) + XcL,n } > 

where Xoptiv) = iulp'^^^ iy)) is the global minimum of 
X^iylp) with respect to p occurring at p°^^{y). The in- 
crement XcL n is ^ fixed number depending on the con- 
fidence level CL and the number of parameters n. It is 
calculated by solving 



E 



CL,n 



Vkix) 



df{x;p) 



dpk 



Denoting A/(a::) = ^Jv{xy E'^^''^ v{x), the global error 
band is approximately given by /max,min(a:) = /(x; p°''^)± 
A/(^) • 



2. Mathematica Code Sample 

A quick way to obtain a statistical er- 
ror band is provided by the following 
Mathematica 5.2 code snippet. The function 
ErrorBand [dat a , model , var s , params , npar , CL] re- 
turns a list containing expressions for /min(a;) and 
/max(2;). The first four arguments data, model, vars 
and params correspond to those in the Mathematica 
function NonlinearRegress: model contains f{x;p), 
vars lists the independent variable(s) x, and params 
specifies the parameters p, optionally together with 
initial values for the search. Here data is a list of tuples 
of the form {xj,yj,Ayj}, i.e. it has an additional column 
specifying the one-standard-deviation errors Ay^ . npar 
is typically set to Length [paramis] so as to generate the 
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global error band, but it must be set to 1 when a local 
error band as discussed in section lA 41 is desired. The 
confidence level CL defaults to 68%, corresponding to 
one standard deviation. If options are given, they are 
passed to NonlinearRegress. 

Please note that ErrorBand [ . . . ] uses an error ma- 



trix which is only a good approximation when the model 
f{x]p) is sufficiently linear in the parameters p and the 
fit is good, cf. eq. (4.62) in ref. 0. A direct calculation of 
the error matrix, as done, for example, by MINUIT, does 
not suffer from this limitation. 



« Statistics ' NonlinearFit ' ; 

« Statistics'NormalDistribution' ; 

ChiSqrThresh[n_, CL_] := 

X /. Solve [CDF [ChiSquareDistribution [n] , x] == CL, x] [[!]]; 

ErrorBand [data_ ,model_ , vars_ ,params_ ,npar_ ,CL_ : . 68269 , options ] : = 

Module [{reg, emat ,grad, thr , err ,p} , 
thr = ChiSqr Thresh [npar , CL] ; 

reg = NonlinearRegress [Drop [data, ■[},{- !}■] ,model , vars ,params , Weights->l/(Last/@data) "2 , 

Regress ionReport->-[BestFit ,BestFitParameters ,EstimatedVariance , AsymptoticCovarianceMatrix} , 
options] ; 

(* extract approximation to error matrix *) 

emat = thr * (AsymptoticCovarianceMatrix/ . reg) [ [1] ] / (EstimatedVariance /. reg); 

(* gradient of fit function with respect to the parameters *) 

p = (Flatten/@(List/(3params)) [[All,l]] ; 

grad = (D [model ,#] &/@p) /. (BestFitParameters/ . reg) ; 

(* calculate parabolic error *) 

err=Sqrt [grad . emat . grad] ; 

■[ (BestFit/ . reg) -err , (BestFit/ . reg) +err} 

]; 



3. Systematic envelope 

Up to here we have neglected uncertainties concerning 
the fixed parameters it, such as gA and While we have 
a good guess about acceptable ranges for their values, we 
have no knowledge about possible correlations of these 
error estimates. Therefore we scan the whole range of 
acceptable values of u on a grid, perform the complete 
statistical error analysis for each choice of u and pick 
the extreme values of /min(a:) and /max (a;). We call the 
resulting region the "systematic envelope" in our plots. 



errors appear on the diagonal of E'^^'^: (Ap^)^ — E'j^^'^ . 
Note that the multi-parameter confidence region is not 
enclosed in the box of single-parameter error bounds 3^ . 

As another important case, consider a local error band, 
formed by plotting the single-parameter confidence inter- 
val for f(x]p) at every x. It is narrower than the corre- 
sponding global band and needs to be interpreted differ- 
ently. Its construction is identical to that of the global 
band, except that now the error threshold is Xcl u ^^'^ 
the error matrix is E'^^'^ . We opt to show such a band in 
fig. llll where we want to read off the error of our nucleon 
mass prediction directly. 



4. Single-parameter errors and local error band 

Confidence intervals for any parameter dependent 
quantity q{p) satisfying P [qmin < q{p°''^) < gmax] = CL 
can be generated in good approximation using the 
method: 

[Qmin, gmax] = { : ivlp) < Xopt (y) + XCL,1 } . 

where now the single-parameter increment Xch i ^'^ 
be employed. For a confidence level CL = 68% corre- 
sponding to a one standard deviation error, Xcl i ~ ^■ 

A simple case of such a parameter dependent quantity 
q{p) is a specific parameter pk itself. With the quadratic 
approximation of from above, the single-parameter 



5. Errors in the pion mass 



The errors Axj in the pion mass can be treated by 
setting up x^ as 



X^iP, Sx) 



N 

E 



f{xj +Sxj;p) -y/ 


2 

-f 


■ dXj ' 




Ay, J 









The minimization and error treatment must now be per- 
formed on n -|- parameters p, 6x. Since we are not 
interested in the joint confidence region of p and 5x, but 
only in the confidence region of p, it makes sense to con- 
tinue using the threshold value Xcl n ^^"^ ^^^P 
minimized with respect to the Sx at all times 0- 
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TABLE V: Two-flavor lattice data for the nucleon mass Mjv taken from IM Holl and selected according to the criteria described 
in appendix |B| The index numbers in the first column comply with those in the compilation 8]. 



no. 


collaboration 


/3 




a [fm] 


L [fm] 


rrin [GeV] 


Mjv [GeV] 


large L 


19 


CP-PACS 


2.1 


0.1382 


0.111 


2.68 


0.5214(21) 


1.2751(82) 


X 


41 


JLQCD 


5.2 


0.1355 


0.098 


1.96 


0.5453(91) 


1.300(23) 


X 


36 


JLQCD 


5.2 


0.1355 


0.099 


1.58 


560fl6~) 


1.412(62) 




31 


JLQCD 


5.2 


0.1355 


0.099 


1.19 


0.655(32) 


1.637(82) 




8 QCDSF 


5.25 


0.13575 


0.092 


2.21 


0.5570(70) 


1.320(20) 


X 


23 


CP-PACS 


2.2 


0.1368 


0.092 


2.22 


0.5946(53) 


1.348(13) 


X 


56 


m 


5.6 


0.1575 


0.085 


2.04 


0.6429(68) 


1.377(19) 


X 


55 


m 


5.6 


0.1575 


0.084 


1.34 


0.660(12) 


1.471(29) 




54 


m 


5.6 


0.1575 


~0.085 


1.19 


0.709(11) 


1.672(38) 




53 


flO] 


5.6 


0.1575 


~0.085 


1.02 


0.832(22) 


1.900(39) 





We can avoid too many parameters by making an ap- 
proximation. Assuming that f{xj + Sxj\p) is approxi- 
mately linear in Sxj for Sxj < Axj and minimizing 
with respect to the Sxj yields Q 

xlsip) = 

y> [fi^j^p) -Vjf 

^tt [Ay.f + [fix, + Ax,/2\p) - fix, - Ax,/2\p)f 

Within our approximation, the same minimum as 

the original, full x^- For hypothesis testing, the number 
of degrees of freedom remains 2N — (TV -I- n) = N — n. 
For confidence regions for p, the x^ increment is Xcl n 
as before. 



APPENDIX B: LATTICE DATA 

The lattice data used for our calculations was taken 
from refs. 0.^3 ^'^'i listed in tabled The conversion 
into physical units was performed setting the Sommer 



scale ro = 0.5 fm Related systematic uncertainties 
need to be estimated by further lattice studies, and are 
not discussed throughout this work. The table is divided 
by horizontal lines into "volume groups" characterized 
by the same value of simulation parameters (3 and k. 
The data have been selected according to the following 
cuts: L > 1 fm, a < 0.15 fm. In addition, the point 
of largest simulation volume in each group must fulfill 
nijr < 0.65 GeV and m^ri > 5. The lattice data from 
[lOj have been generated with a standard, unimproved 
Wilson fermion action. Even though, lattice artefacts 
are shown ^3 to be small for the points we select here, 
in virtue of the fine lattice spacing. For points 53 and 54, 
no calculation has been performed for the Sommer radius 
To/a, which is needed to determine the lattice spacing 
a. In accordance with ref. |lOj [. a for these two points is 
copied from point 56. 
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